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O ■ ABSTRACT 
O 
(N 

Q ■ This is the first attempt to construct detailed X-ray spectra of clusters of 

galaxies from the results of high-resolution hydrodynamic simulations and sim- 
ulate X-ray observations in order to study velocity fields of the intracluster 
medium (ICM). The hydrodynamic simulations are based on the recently pro- 
posed tsunami model, in which cluster cores are affected by bulk motions of the 
> ■ ICM and turbulence is produced. We note that most other solutions of the cooling 

flow problem also involve the generation of turbulence in cluster cores. From the 
' mock X-ray observations with Astro-E2 XRS, we find that turbulent motion of 

^ the ICM in cluster cores could be detected with the satellite. The Doppler shifts 

^ ■ of the metal lines could be used to discriminate among turbulence models. The 
O ■ 

' gas velocities measured through the mock observations are consistent with the 

Q^' line-emission weighted values inferred directly from hydrodynamic simulations. 

O 

^ I Subject headings: cooling flows — galaxies: clusters: general — turbulence — X- 



in 



rays: galaxies: clusters 



H , 1. Introduction 

The advent of Astro-ES (Inoue 2003) will enable us to directly measure velocity fields of 
the intracluster medium (ICM) in clusters of galaxies for the first time. The calorimeter (X- 
Ray Spectrometer; XRS) of Astro- E2 has an excellent spectroscopic resolving power (Kelley 
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2004; Cottam et al. 2004), and it could detect the bulk gas motion by observing the energy 
shift of metal lines and the turbulence by observing broadened metal lines. On the other 
hand, many hydrodynamic simulations have been performed to study the motion of the 
ICM (e.g. Evrard 1990; Takizawa 1999; Yoshikawa et al. 2000; Ricker & Sarazin 2001). With 
Astro-E2, we could confirm the results of those hydrodynamic simulations. However, the 
direct comparison of X-ray observations to hydrodynamic simulations is not always simple, 
because of complex such as instrumental responses. 

In this paper, we construct the X-ray spectra from the results of hydrodynamic simu- 
lations of cluster cores affected by 'tsunamis' (Fujita, Matsumoto, & Wada 2004, hereafter 
Paper I, see also Fujita, Suzuki, & Wada 2004) and we simulate observations with Astro- 
E2 XRS. In the tsunami model, large-scale bulk gas motions in the ICM, which are called 
tsunamis, induce fully-developed turbulence in cluster cores. The ultimate spectroscopic ca- 
pabihty of XRS with an energy resolution of ~ 6 eV at 6 keV for extended sources as well as 
point sources could give us a new approach to detect the turbulence in the ICM. Because of 
the kinetic energy of turbulence and the hot gas brought into the core by turbulent mixing, 
the radiative heating of the core is suppressed. Thus, this model could solve the so-called 
'cooling flow problem' (Makishima et al. 2001; Peterson et al. 2001; Kaastra et al. 2001; 
Tamura et al. 2001). 

From the mock observations with Astro-E2, we measure velocity flelds of the ICM and 
compare them with the emission weighted values inferred directly from the hydrodynamic 
simulations. We investigate whether the former is consistent with the latter. Most observa- 
tions of clusters scheduled as Astro-E2 performance verification targets in the first 6 months 
are planned pointing at their centers^; one of the reasons is that the effective area of the 
XRS is relatively small and objects must be bright enough to be observed. Thus, it would be 
useful to simulate the observations focused on cluster cores. It should be noted that other 
solutions of the coohng flow problem besides the tsunamis, especially those based on AGN 
activities, also predict similar level of turbulence in cluster cores, which some observations 
have already suggested (e.g, Blanton et al. 2001; Churazov et al. 2002; Briiggen & Kaiser 
2002; Kim & Narayan 2003; Kaiser & Binney 2003; Fabian et al. 2003; Soker et al. 2004). 
In this paper, we assume that VLq = 0.27, A = 0.73, and Hq = 70 km s~^ Mpc~^. 
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2. Construction of X-Ray Spectra 

We use the results of hydrodynamic simulations done in Paper I. These simulations are 
performed using a two-dimensional (cylindrically symmetric) nested grid code (Matsumoto 
& Hanawa 2003), and the coordinates are represented by {R,z). Since the structure of 
turbulence in the axi-symmetric coordinate might be different from that in the fully three- 
dimensional coordinate, it is ideal to perform three-dimensional calculations. However, a 
high spatial resolution with a large dynamic range is also crucial to simulate a turbulent 
medium, especially for the tsunami model. Therefore we here restrict the calculations to 
two dimensions. The cluster center corresponds to {R,z) = (0,0). Although the maximum 
resolution of the simulations is achieved on the level of grids of 22 pc, we use the simulation 
results on the level of the grids of 1.4 kpc, because the turbulent velocity is generally smaller 
on smaller scales and it helps us reduce the number of computational grid points for which 
we calculate the spectra. We will argue that the effect of the lower resolution on the X-ray 
spectra can be ignored in §4. Note that in the nested grid code we used, calculations are 
performed on all grid levels simultaneously. At the cluster center, for example, the solutions 
(density, temperature, etc.) are obtained on the seven different resolutions. The results of 
the lower resolution grids are just the coarsened ones of the higher resolution grids. 

We use the simulation results for 256 x 512 computational grid points. For each grid 
point, we calculate the X-ray spectrum using a single bapec thermal model and the data 
simulation command fakeit in the XSPEC package (version 11.3.1). The bapec model 
is the same as the apec model but includes a parameter for turbulent velocity broadening 
to emission lines. The input parameters of the spectrum for the i-th grid point are the 
temperature (Tj), the metal abundance (Zi), the gas velocity along the line of sight (fj), and 
the normalization. The turbulent velocity for each grid point is assumed to be zero, that is, 
the gas velocity is uniform within a single grid point. We used the Astro-E2 XRS response 
file (xrs_aol .rmf) and the auxiliary file (xrs_onaxis_open_aol .arf) provided by the Astro- 
E2 team as planning tools for AO-1 proposals^. Weighting the three-dimensional volumes 
of the grid points, the spectra are summed up by the FTOOLS manipulation task mathpha. 
In order to avoid the error caused by small photon counts, we multiply the photon counts 
for each grid point by 100000. After the summation by mathpha, the total photon counts 
are divided by 100000 using the f calc task in the FTOOLS. Since f calc discard fractions 
in the calculation of the photon counts, we add 50000 to the counts before the division in 
order to round off. For the resultant spectrum, the error of photon counts in each energy 
bin is assumed to be the square root of the photon counts in the bin. Background emission 
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and Galactic absorption are ignored. We do not consider the effect of resonance scattering 
for metal emission lines, because the calculation of multi-dimensional radiative transfer is 
required and it is beyond the scope of this paper. 

3. Results 

In paper I, we calculated the evolution of the cool core of a cluster for R < 300 kpc and 
1^1 < 300 kpc. We approximated the bulk gas motions in a cluster by plane wave-like velocity 
perturbations in the ^-direction represented by Svz = acg sm(2n Cst/X) at 2; = —345 kpc, 
where Cg is the initial sound velocity, and A is the wave length. The waves propagate in the 
2;-direction. At t = 0, the cluster is isothermal (T = 7 keV). Metal abundance is uniform 
and Z = 0.3 Zq. Wc construct X-ray spectra for two models with different A (Table 1). 
The model of a = 0.3 and A = 100 kpc (Model A) was studied in Paper I. The model of 
a — 0.3 and A = 1000 kpc (Model B) is newly studied and the detailed analysis of the results 
will be discussed elsewhere (Matsumoto, Fujita, & Wada 2004, in preparation). Compared 
to Model A, Rayleigh-Taylor and Kelvin- Helmholtz instabihties develop on larger scales in 
Model B. This makes turbulent mixing and heating more effective. 

We assume that the model cluster is at redshifts of 0.01, 0.04, and 0.08, although the 
temperature outside of the core is ~ 7 keV and there is no such a high-temperature cluster 
at redshift of 0.01. Since the field of view of Astro-E2 XRS is 2.9' x 2.9', we consider the 
X-ray emission within a radius of 1.5' from the cluster center. At redshifts of 0.01, 0.04, 
and 0.08, the angle of 1.5' corresponds to 18.5, 71.3, and 136 kpc, respectively. We 'observe' 
the model cluster along the 2;-axis, that is, the line of sight is assumed to be parallel to the 
z-axis. This is because the waves were injected along that axis, and thus the ICM velocity in 
the 2;-direction is much larger than that in the i?-direction. Thus, we sum up the spectra of 
the ICM in individual grid points for R < 18.5, 71.3, and 136 kpc and \z\ < 300 kpc for the 
cluster at redshifts of 0.01, 0.04, and 0.08, respectively. This means that for a cluster at a 
larger redshift, we observe a larger area. The assumed exposure time is 50 ks. Using XSPEC, 
we fit the summed spectra with a single bapec model. The spectra were grouped to have 
a minimum of 50 counts per bin. We limited the energy range to 5-10 keV (at the cluster- 
rest frames) that includes Fc-K lines (~ 6.7 keV). Free parameters in the bapec model are 
the temperature (T), the metal abundance (Z), the average velocity (f^.sp)) the turbulent 
velocity ((T^ ,5^), and the normahzation. Wc choose t = 3.3 Gyr for Model A; the temperature 
distril)ution at that time is shown in Paper I. We choose t = 5.0 Gyr for Model B.'"^ For 
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Model A, the gas temperature in at least one of the grid points reaches zero at t — 3.3 Gyr 
(Paper I), while for Model B, it reaches zero at t = 6.2 Gyr. Note that the time-scale of 
6.2 Gyr is comparable to the typical age of clusters (e.g. Kitayama & Suto 1996), and thus 
radiative cooling is almost completely suppressed in Model B. 

Figure 1 shows the simulated X-ray spectrum around 6.7 keV Fe-K hues for Model B at 
redshift of 0.04 and the result of the fit. For comparison, the spectrum when Vz,sp = and 
o'z,sp = is also shown. As can be seen, the two spectra (i.e. with and without turbulent 
motion) are remarkably different. The results of the spectral fits for all models are shown 
in Table 2. Errors on fitted spectral parameters are given at the 90% confidence level. 
In Figure 2, we show Vz,sp and az,sp- Most models require non-zero turbulent velocity. The 
values of are generally small and the fits are good. The results show that the fine shift and 
turbulence will be detected with Astro-E2 XRS if Vz,sp ^ lOOkms"^ and a^^sp ^ lOOkms"^, 
because Vz,sp — and az,sp — are rejected, respectively (Table 2). However, it should 
be noted that our results correspond to an ideal case. For example, we did not consider 
possible systematic errors remaining even after calibrations. Moreover, if the line of sight 
is perpendicular to the ;z-axis, observed average and turbulent velocities are smaller. Thus, 
turbulence could not be detected for all clusters. 

For comparison, we present the luminosity- weighted turbulent velocities in the z-direction 
that are directly derived from the results of hydrodynamic simulations. They are given by 



<^zA = 



'^i{Vz,i-Vz,iyLi 



where Vz^i is the gas velocity in the ^-direction, and Lj is the luminosity of the i-th grid point. 
As for the luminosity, we consider both the bolometric luminosity, Lboi.i, and the Fe-K-line 
luminosity Lpci- Using XSPEC and the bapec model, we found that the emissivity to be 
used for the latter is approximately given by 

epe = 5.559 X 10"^VenH^ ergs cm'^s"^ (2) 

X ^Mi^^:^^lI±Il!- + l-4.215T + 0.2492T2- 0.005255 T3 , 
13.14+ (0.00169 + T)2 J ' 

for the 6.6-6.72 keV band, where T is the ICM temperature (keV), Z is the metal abundance 
{Zq), and Tif. and tih are the electron density and the hydrogen density (cm^'^), respectively. 
This approximation holds within a few percent for 1.5 < T < 15 keV (even at T ~ 1 keV, 
the error is ^ 10 %). We ignore the Fe-K emission from the ICM of T < 1.5 keV, because 
epe ~ and the mass of such gas is small in our calculations. In equation (1), the average 
gas velocity in the z-direction, Vz^u is given by 

_ Y^i Vz,iLi , , 
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We refer to az,i and Vz,i weighted by the bolometric luminosities as az,i,boi and Vz,i,boi, respec- 
tively, and those weighted by the Fc-K luminosities as (J^,i,Fe and Vz^i^fc, respectively. The 
velocities are shown in Figure 2 and az,i,Fe ^-nd Vz,i,Fe are also shown in Table 2. 

4. Discussion 

In Figure 1, the line profiles appear nearly symmetric. This is because the probability 
distribution function for Vz is almost symmetric (Matsumoto et al. 2004, in preparation). 
Recently, Inogamov & Sunyaev (2003) investigated line emission from turbulent gas in detail 
and indicated that line profiles could be very complicated; asymmetry would be observed in 
the X-ray spectra. However, for the turbulence we considered, the expected asymmetry is 
too small to be detected with the energy resolution of Astro-E2 XRS. On the other hand, the 
asymmetric injection of tsunamis is observed as a large offset of Vz ~ —140 km s^^ (Model B; 
Table 2), because the asymmetric tsunamis (e.g. those created by cluster mergers) sometimes 
induce not only the turbulence but also the bulk motion (or oscillation) of the cool core; the 
latter is unhkely to be induced by symmetric jets accompanied by AGN activities. Therefore, 
the relative motion between the cD galaxy (stars) and the strongly turbulent core (gas) in a 
cluster could be a clue to confirm the tsunami model as well as the detection of turbulence 
in clusters without AGN activities. 

Figure 2 shows that Vz,sp and Vz,i,Fe are consistent, and also show that (7z,sp and (7z,i,Fe are 
consistent. These mean that velocity fields obtained through hydrodynamic simulations can 
directly be compared with X-ray observations without resorting to mock observations like the 
one we did here, if the velocities are weighted by line-emission. On the other hand, near the 
cluster center, which means at smaller redshift, in Model A, Vz^sp and c^z^sp are not consistent 
with Vz,i,boi and az,i,boi, respectively (Figure 2). However, they are consistent in Model B. 
The difference between Model A and B is because radiative coohng proceeds more and the 
gas temperature at the cluster center is lower for Model A ai t — 3.3 Gyr than for Model B 
a,t t — 5.0 Gyr. In Model A, the temperature near the cluster center is < 3 keV, which is 
much smaller than the temperature outside of the core (~ 7 keV) and the Fe-K emission is 
relatively weak. Since cooler gas tends to have larger bulk and turbulent velocities, the gases 
with larger velocities are less weighted by the Fe-K luminosities. Therefore, in Model A, the 
velocities weighted by the bolometric luminosities (|f2,/,6o«| and (Tz,i,boi) larger than those 
weighted by the Fc-K luminosities (It'^./.Fcl and a^j pe) near the cluster center. On the other 
hand, in model B, the gas temperature in the central region of the cluster is ~ 4 keV and 
the temperature gradient near the center is small compared to that in Model A. Therefore, 
Vz,i,boi {(7z,i,boi) and Vz,i,Fe {(7z,i,Fe) are almost the same. 
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In Figure 2, the luminosity- weighted turbulent velocities, az,i,boi and (Jz^i^Fe, increase 
toward the cluster center (toward smaller redshift). This means that the turbulence produced 
by tsunamis is more detectable at the cluster center. 

As mentioned in §2, we used a low resolution grid. However, we expect that the effect of 
the coarsening on the Vz,sp and az,sp is small, because Vz,i,Fe and az,i,Fe are not much different 
in between the low resolution grid and the highest one; the difference is much smaller than 

the errors of v^^sp and (Tz,sp presented in Table 2 and Figure 2. Since Vz^sp and Vz^ifc {<^z,sp 
and (T^,i,Fc) are almost same as discussed above, Vz^sp {<^z,sp) will not change much even for 
the finest grid. 

We also 'observed' the cluster from the direction perpendicular to the ^-axis and we call 
this direction x. Because of the symmetry we assumed, the average velocities, v^., are zero. 
The turbulent velocities derived from the mock observations are (Jx^sp ~ 60 and 120 km s~^ 
for Model A and B, respectively. Consistency between ax,sp and those weighted by Fe-K 
lines, ax,i,Fe, is also good. 

5. Conclusions 

We have constructed X-ray spectra from the results of hydrodynamic simulations of 
clusters of galaxies based on the tsunami model, in which turbulence is created in cluster 
cores. Similar levels of turbulence are also predicted by many of other heating models of 
cool cores (e.g. motion of bubbles created by AGN activities; Churazov et al. 2002; Briiggen 
& Kaiser 2002). In particular, we focus on the effect of velocity fields of the ICM on the 
spectra. We simulate X-ray observations with Astro- E2 XRS and find that velocity fields 
in cluster cores could be revealed with the sateUite. The motion of the cool core could 
be used to discriminate among turbulence models We show that the gas velocities derived 
through the mock observations are consistent with the Fe-K line emission weighted values 
inferred directly from hydrodynamic simulations. The technique developed here could easily 
be applied to the comparison between results of various hydrodynamic simulations and those 
of near-future observations with Astro-E2 and others {Constellation- X, XEUS, NeXT). 
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Table 1. Model Parameters 



Model 


a 


X (kpc) 


* (Gyr) 


A 


0.3 


100 


3.3 


B 


0.3 


1000 


5.0 



- 11 - 



Table 2. Fitting Results 



Model 
(Redshift) 


T 
(kcV) 


Z 

(Zo) 


^2 , sp 

(km s^i) 


(^z .sp 

(km s^i) 




(km s^i) 


(km s^i) 


(5-10 kcV)'= 
( erg cm^2 s^i) 


A (0.01) 
A (0.04) 

A (0.08) 
B (0.01) 
B (0.04) 
B (0.08) 


2-80_o.o5 

4 r;p;+0.20 
*-OO_0.19 

5 71+0.41 

' -"--0.47 

4 04+"-" 

, 01+0.36 
'*-°^-0.33 


0.26lHi 

0.28t°:0^2 

sn+0 03 

2q+"-"i 

'^••^'^-0.02 
29+° °^ 


10+15 

9+13 
^-16 

t; + 25 
"^-28 

-147l?i 

-123l2? 

— 1 

^'^'-'-36 


256ti? 

230t2? 
20311 


217.6/679 
122.8/241 

170.4/95 
163.8/724 
128.5/283 

88.2/104 


12 
-2 

2 

-150 
-128 
-158 


124 
83 

80 
237 
214 
201 


1.1 X 10-1° 

3.2 X 10-" 
1.2 X 10-" 
1.1 X 10-1° 
3.8 X 10-11 
1.4 X 10-11 



''Average velocity weighted by the Fe-K lines luminosities 
''Turbulent velocity weighted by the Fe-K lines luminosities 
"^X-ray flux in the 5-10 keV band 
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T ' 1 ' 1 ■ r 




channel energy (keV) 

Fig. 1. — Upper panel shows the simulated X-ray spectrum observed with Astro-E2 XRS for 
Model B at redshift of 0.04 (crosses) and the best-fit model (solid line). The spectrum when 
Vz,sp = and a^^sp = is also shown (dashed line). The spectra are shown at the cluster- rest 
frame. The lower panel plots the residuals divided by the la errors. 
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Fig. 2. — (a) Average gas velocities for different redsfiifts. Radii corresponding to the field 
of view of XRS (R) are also shown. Dotted lines indicatev z^i,boi andsolid lines indicate Vz,i,Fe- 
Filled circles indicate Vz,sp- Bold hnes and large circles correspond to Model A and thin 
lines and small circles correspond to Model B. (b) Turbulent velocities for different redshifts. 
Dotted lines indicate (Tz,i,boi and solid lines indicate (Tz,i,Fc- Filled circles indicate (Xz^sp- Bold 
lines and large circles correspond to Model A and thin lines and small circles correspond to 
Model B. 



